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Abstract 

This  report  contains  a  first  order  analysis  of  the 
supercritical  flow  in  an  overfall  based  on  the  use  of  the 
stream  line  function  as  the  basic  dependent  variable.   It 
presents  approximate  equations  for  the  nappes;  and  the 
formula 


f*(i;e) 


1 


1 


where 


/IT^   >-'l+£2(.772) 


ga 
u 


g  is  the  gravitational  acceleration,  a  is  the  upstream  depth 
and  u  is  the  velocity  there.   This  formula  is  a  good  approxi- 
mation to  the  brink  depth  ratio  b/a  for  the  range  0  <  e  <  1. 


Ill 


1.  Introduction 

This  report  is  concerned  with  the  supercritical  flow  in  a 
two-dimensional  overfall.   The  analysis  is  based  on  the  use  of  the 
streeun  line  function 

y  =  f(x,7)  . 

This  function  is  such  that  if  the  stream  function  is  ^(x,y)  so  that 
the  stream  lines  are  given  implicitly  by 

^(x,y)  =  7  , 
then 

^{x,f(x,7))  37 

is  identically  satisfied. 

if  the  acceleration  due  to  gravity  is  g,  if  the  upstream  depth 
is  a,  and  if  the  velocity  there  is  u,  then  the  flow  is  said  to  be 
critical  if  e  =  ga/u  is  equal  to  1,   We  suppose  that  e  is  small  so 
that  the  upstream  velocity  (u  >>  yga)  of  the  flow  is  supercritical. 
Such  a  flow  can  prevail  in  a  region  of  space  where  g  is  small  compared 
with  a  and  u.   It  can  obviously  also  prevail  if  u  is  sufficiently 
large,  or  if  a  is  so  small  that  the  overfall  is  nearly  a  falling 
sheet. 

We  show  below  that  the  stream  line  function  must  satisfy  a 
nonlinear  partial  differential  equation  defined  in  an  infinite  strip, 
and  that  it  must  satisfy  nonlinear  boundary  conditions  along  the 
side's  of  the  strip.   Since  these  boundary  conditions  involve  the 


small  parameter  e,  we  assume  that  T{x,-y)    can  be  approximated,  at 
least  asymptotically,  by  an  expansion  in  integral  powers  of  e. 

The  development  which  follows  is  confined  to  an  analysis  and 
discussion  of  the  first  order  approximation  and  what  it  yields  with 
respect  to  the  shapes  of  the  free  surfaces  and  the  brink  depth  of 
the  overfall.   A  formula  which  is  derived  for  the  brink  depth  ratio, 
namely. 


f*(i;£)  =  ^ 


/T+2l 


+ 


'l  +  E^(.3l4) 


l  +  e^(.772) 


appears  to  be  a  good  approximation  in  the  range  0  _^  e  ^  !•   According 
to  this  formula  the  brink  depth  ratio  for  critical  flow  is  .719* 
This  is  close  to  .708,  the  average  of  six  approximations  based  on 
different  procedures  developed  by  various  investigators. 


2.  Formulation 

A  two-dimensional  waterfall  or  free  overfall  is  defined  by- 
its  flow  in  a  plane  perpendicular  to  its  linear  edge.   Figure  2.1 


OD 


D 


\ 


►  X 


Fig.  2.1 


shows  such  a  flow  referred  to  a  rectangular  coordinate  system  whose 
origin  is  at  the  brink.  The  negative  part  of  the  x-axis  coincides 
with  the  sill  or  bottom  L-,  and  the  y-axis  is  positive  upwards.  The 
domain  D  of  the  cross  sectional  flow  is  bounded  by  L-,;  and  the  curves 
r-i  and  Tp  which  define  the  free  surfaces  or  nappes  of  the  overfall. 
The  difficult  basic  problem  is  to  deduce  the  -shapes  of  Tn  aJ^d 
To  a-s  consequences  of  the  assxamptions  that  gravitational  attraction 
is  the  only  body  force  acting  and  that  the  flow  is  a  steady  irrota- 
tional  one  of  an  inviscid,  incompressible  fluid  whose  density  5  is 


constajit.   If  _v  is  the  velocity  vector,  conservation  of  mass  requires 
V'V  =  0  and  this  with  the  vorticity  condition,  Vx  v  =  0,  insures  the 
existence  of  a  function  of  the  cojnplex  variable  z  =  x+  iy,  namely  the 
complex  velocity  potential, 

(2.1)  w  =  c|)(x,y)  +i^(x,y)  =  F(z) 

which  is  analytic  for  z  in  D.   This  function  is  such  that  the  gradient 
of  its  real  part  gives  the  velocity  vector  v  =  V  Re  F(z)  =  Vcf);  and 
such  that  its  imaginary  part,  i;i/(x,y)  =  Im  F(z),  is  the  stream  function 
for  the  flow.   In  other  words,  the  horizontal  and  vertical  velocity 
components  are  given  respectively  by  the  real  and  Imaginary  parts  of 


F'(z)  =  c|)-i^   =?//+i(j) 
^  ^    ^x    ^x    ^y    ^y 

If  n  is  the  unit  outward  normal  to  the  boundary  of  D  and  if  s 
denotes  arc  length  along  the  boundary,  then  the  kinematic  boundary 
condition  requires  that  the  normal  velocity 

*n  =  ^s 

must  be  zero  along  L-, ,  T-,  and  Pp.   It  follows  from  this,  by  integra- 
tion with  respect  to  s,  that  L-,  plus  r-,  constitutes  a  stream  line 
which  we  can  assume  is  defined  by  the  level  curve 

^(x,y)  =  0  . 


With  this,  the  upper  nappe  curve  To  is  defined  by  level  curve 


^(x,y)  =  h  =  au 

where  a  is  the  upstream  depth  of  the  flow  at  x  =  -co,  and  u  is  the 
velocity  there.   In  addition  to  the  kinematic  boiindary  condition, 
the  energy  integral  of  the  momentum  equation 

6.y_  /  ^Z       "N 

^  dt  =  5(^^+  v.vvj=  5V.VV 

=  -p  [Vv«v  -V  XVX  v] 


_  vv. V  =  -Vp  -  g5Vy 


namely. 


(2.2)  ^^+<^^+^+2gy  =   c 

must  he  satisfied  along  r-,  and  Tp.   In  this  equation  g  denotes  the 
gravitational  acceleration  and  p  denotes  the  pressure  in  the  fluid. 
Since  the  flow  is  supposed  to  be  steady  the  constant  c  remains  the 
same  no  matter  where  the  left-hand  side  is  evaluated.   Therefore  with 

<l)^(-oo.y)  =  u 

(f>y(-oo.y)  =  0 


p(-oo,a)  =  p^ 


we  have 


(2.3)  c  =  u^+2ga+-^ 


In  what  follows  we  assume  that  the  pressure  on  each  nappe  is  constsuit 
and  equal  to  p^. 


Although  the  domain  D  is  not  known  at  the  outset  we  can  work 
with  a  fixed  domain  in  the  w  =  '^  +  lip   plane  if  we  introduce  the  inverse 
function 


(2.4) 


z  =  x((|),^)+  iy(({),^)  =  G(w)  . 


Under  the  mapping  from  the  z-plane  to  the  w-plane,  D  is  mapped  into 
the  strip  defined  by 


0<-^  =  Imw<h 


au 


-co    <({)   =  Rew<oo    , 


and  we  can  choose  the  origin  of  the  w-plane  to  be  the  image  of  the 
origin  of  the  z-plane.   With  the  origin  fixed  in  this  way  the  images 
of  L-, ,  r-i  and  Po  are  respectively  M-, ,  Z-,  and  2p  as  shown  in  Figure 
2.2.   Since  there  are  no  stagnation  points  in  the  overfall,  the 


Fig.  2.2 


function  G(w)  must  be  analytic  in  S  and  we  have 


§  =  G'(w)  =  x^+iy^  =  -^^^-^y^ 

dz  "^    ^x    ^y 


from  which 


and 


G'(w)G'(w)  =  |G'(w)|2  =  x?  +  y?  =   o^  o 

'Px  +  'f  y 

.   ^  Re  G ' ( w )      I   ^  Im  G' (w) 
^^  "  |G'(w)r^  '    ^y  ~    |G'(w)|'^  * 

These  relations  with  (2.2)  and  (2.3)  show  that  G(w)  must  satisfy  the 
following  boundary  conditions: 

Im  G'  ((f))  =  0  ,  ^    <   0      ; 

^  +  2g  Im  G((|))  =  -u^  +  2ga  ,  (f)  >  0   ; 

|G'(4)r 

+  2g  Im  G((l)+ih)  =  u  +  2ga  ,        -oo  <  ({)  <  oo  . 


|G'('l'+ih)l 
We  can  now  use  the  transformation 

w  =  m((^ ) 

to  map  S  into  any  domain  we  think  may  be  helpful  for  the  analysis. 
(For  a  strictly  numerical  analysis,  a  mapping  of  S  into  a  rectangle 
offers  some  advantages.)   In  a  previous  report,  [1],  the  author  used 

w  =  m(C)  =  ^  ^n  (-  ^)  . 

This  maps  S  into  the  upper  half  of  the  ^-plane  and  gives 


G(w)  =  G[|  in    (-  ^)]  =  o(tJ  =  ai+  MO    • 

In  [1]  a  nonlinear  integral  equation  was  derived  for  the  determina- 
tion of  ou(C)  and  then  approximations  were  found  for  G(w). 

In  the  beginning  the  problem  was  to  determine  a  function  F(z) 
analytic  in  an  unknown  domain  D  but  subject  to  two  boundary  condi- 
tions along  the  \inknown  part  of  the  boundary  of  D.   We  have  seen  how 
the  dual  boundary  conditions  enable  us  to  transform  the  problem  into 
one  for  the  determination  of  a  function  o(0  analytic  in  a  known 
domain  but  subject  to  a  nonlinear  boundary  condition.   The  nature  of 
the  nonlinearlty  is  such  that  the  probability  of  finding  o(C)  in 
closed  form  is  practically  nil.   We  can  only  expect  to  find  approxima- 
tions for  o(C)  and  hence  G(w).   Note  however  that  such  approximations 
present  approximations  for  the  equations  of  r-j^  and  Tp  ih  terms  of 
the  parameter  (j).   For  example,  if  G(w)  is  an  approximation  to  G(w) 
then  the  approximate  parametric  equations  for  y^   are 

X  =  Re  S((|)+  ih) 
y  =  Im  G(<j)+  ih)  . 

This  parametric  representation  is  not  an  advantage  when  it  comes  to 
calculating  the  depth  of  the  liquid  at  some  prescribed  distance  from 
the  edge  of  the  overfall.   One  such  depth,  the  brink  depth  b  in 
Figure  2.1,  is  perhaps  the  most  important  single  numerical  value 
which  is  desired  from  an  analytis  of  the  overfall.   In  order  to  calcu- 
late this  depth  we  must  at  least  approximate  the  solution  of  the 
trajiscendental  equation 


8 


0  =  Re  Sf(({)+  ih)  , 

This  kind  of  burden  is  of  course  inherent  in  the  above  kind  of 
analysis  and  it  must  be  carried  whenever  we  wish  to  estimate  hydro- 
dynamic  q-Ucintities  at  prescribed  points  in  the  overfall. 

Instead  of  using  cj)  and  tf/   as  independent  variables  for  the  flow 
we  can  introduce  the  dimensionless  quantities  ^  =  x/a,  f  =  y/a, 
^(x,y)/au  =  T;i/(a|,  af  )/au  =  -q;  and  then  use  |  and  t\   as  independent 
variables  with  f  as  the  basic  dependent  quantity.   The  equation 

(2.5)  ^(a^,  af )  =  Tjau 

defines  the  stream  lines  implicitly.   However,  since  there  is  no 
stagnation  point  in  the  flow  we  can  solve  (2.5)  for  af.   This  gives 

(2.6)  |=f(4,n) 

which,  since  it  presents  the  stream  lines  in  explicit  form,  can  be 
called  the  stream  line  function.   The  range  of  |  is  from  minus 
infinity  to  plus  infinity  and  the  range  of  t]  is  from  zero  to  one. 
The  dimensionless  velocity  components 

^1     ^y(^^y) 

~u       u      ^'^y 
and 


Vg    i^^{^,y) 


u  u  'X 


-ar]. 


can  be  found  in  terms  of  f  by  differentiating  (2.6)  with  respect  to 
y  and  x.   We  have 


T   =  f  (^.ii)t1. 


(2.7) 


n 


which  give  the  dimensionless  horizontal  component  of  velocity 

(2.8)  —L=ar]      =  -^r-r4 r 

^    ^  u     'y   f^(?,TiJ 

and  the  dimensionless  vertical  component 

The  equation  for  f  can  be  foiind  by  recalling  that 


'xx   'yy 


It  turns  out  that  f  must  satisfy  the  nonlinear  equation 
(2.10)  4f£l]   1|^(^)=  0 


1  a  / 

^1+  f|" 

2  ^\ 

or  what  is  the  saime  equation 


(2.11)  (1+  f^)f   +  f^f..  -  2f.f  f.   =0, 


in  the  strip  S*  shown  in  Figure  2.3.   Along  L-,,  f{^,r])   must  satisfy 
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I 


s* 


* 

^2 


^^- 


(2.12) 


Fig.  2.3 
f(e,0)  -  0 


Along  r-|  and  Toj;  f(^>Tl)  niust  satisfy  the  condition  which  comes  from 
(2.2),  namely. 


(2.13) 
where 


1+  ff +  (2ef  -  1  -  26 )f^ 


e  -  -^ 
u 


=  0 


is  small  when  the  flow  is  supercritical. 

If  e  is  zero  the  equation  (2.11)  and  the  boundary  conditions 
(2.12)  and  (2.13)  are  satisfied  by 


11 


This  corresponds  to  a  simple  parallel  flow  with  velocity  components 


V-,  =  u 


and  stream  lines 


v^  =  0 


Yo  =  ari 


If  e  is  small  a  plausible  assumption  is  that  f(^,ri)  =   f(4,ri;e)  can  be 
expanded  in  the  form 


(2.14)       y^^-J^^^^  =  f(e,Ti;e)  =  Ti  +  XZ:  e^d.n) 


which  is  at  least  asymptotically  valid  provided  ^  is  bounded  away  from 
plus  infinity  where  f  is  infinite.   Our  immediate  objective  now  is  to 
analyze  some  of  the  consequences  of  the  above  ass^umption. 


12 


3 .  First  Order  Approximation 

If  we  substitute  (2.l4)  in  (2.11)  and  the  boundary  conditions 
(2.12),  (2.13)  aind  then  equate  coefficients  of  like  powers  of  e  we 
are  led  to  a  sequence  of  linear  elliptic  boundary  value  problems  for 
the  fj^'s. 


(3.1) 


The  function  f-,  must  satisfy  the  equation 


'i^I^^'^^^St  =   ° 


-a:i  <  ^  <  OD , 
0  <  Ti  <  1 


and  the  boundary  conditions 


(3.2) 


r  f-^(^,o)  =  0  , 


i  <  0  , 

0  <  e  , 

-GO   <  ^  <  00 


The  function  fp  must  satisfy  the  equation 


(3.3) 


f2^^(I.Ti)  +  f2^^  =  2(f^^f^^^+f^^f^^^)  , 


■OO  <  ^  -^03, 
0  <  T]  <  1 


and  the  boundary  conditions 


(3.4) 


f2,(^.0) 


=  0  , 


1  ^2 


f^{i,0)+±   f^^(^,0)+| 


^271(^,1)  =  fi(e,i)  +1  fi^(e,i) 


^  <  0  , 
e  >  0  , 

-OO  <  ^  <  OO 
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If  we  introduce  the  complex  variable 

C  =  ^  +  iTi  , 
then 

f-L(^,Ti)  =  Re  F^{^) 

where  F-j^(0  is  analytic  when  C,   is  in  S* .   Along  the  boundary  of  S* 
F-,(0  must  satisfy  the  conditions 

Re  f[u)  -  0  ;  e  <  0  , 

Im  f|(?)  =  1  ,  ^  >  0  , 

Im  F-,(|+i)  =  0  ,  -oo  <  ^  <  oo  . 

Since  the  velocity  at  the  brink  of  the  overfall  is  neither  zero  nor 
infinite,  F-.{C,)   must  be  bounded  at  (^  =  0.   We  suppose  that  the  singu- 
lar behavior  of  F-,(C)  as  ^  -» +oo  is  as  mild  as  possible.   An  analysis 
then  shows  that  the  function  which  satisfies  the  above  conditions  is 


(3-5)  FAO'il      <n[    '/l-e'^^-  1    IdA    ■ 


0  \V  1 


-e-^^+1 


The  first  order  approximation  to  the  stream  line  function  is 


y-,  (^,-n;e ) 

(3.6)       ^—^ =  ^+Ef^{^,r]} 


v/l-e-^^-  1 


=  Ti+:^Re/   ^n         ±   )  h^ 

0     L'l-e-^+l. 
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As  far  as  we  know  the  integral  in  (3.6)  is  neither  tabulated  nor 
expressible  as  a  finite  combination  of  tabulated  functions.   However, 
various  pertinent  expansions  and  estimations  can  be  found  without 
much  trouble. 

The  representation  (3*6)  can  be  expressed  in  the  form 


y 

"a 


A 


1=  ^  +  l:Iief  [in    (-e-^^)-2   ^n  (Jl-e-^^+1)] 
0 

/    r-Tr(A-i)   -2    (^n   2]dA 
0 

oo        ^  I ^ 


dA 


=  n  + 


^  \ 


dA    )    . 


0 


<0D 


Since 


+  s  ;     <n  I  ^1-^1''^+^  ,  d* 


-    /".n  ^iHi^)  ..  =   -  1/     ,n  (/E^)  « 


2     r     ^n    (1-g^) 
TT  J  1+a 

0 


da 


=  -  i  ^n22  +  I 
TT  6 


we   have 


ttC 


+    (ttI-  2    ^n   2JC 


—  =  Ti  +-^  Re 
a  '     TT 


-|.n^..g../.n(£i^) 


. 


dA 


y 


For  Re  C  =  ^  positive  and  large,  the  Integral  in  the  last  expressi 


on 
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is    small.      Therefore,    far   downstream,    we   have   the   approximation 


sa 
y-^  ^  a.T]  +  ^^ 

U    TT 


r  IT       /.2      ,     IT  2  I  ^       „  r^\, 

-    n   i     +-^  T\    -   [2   Zn  2)1 


-  TT^   -  I  in22  +  I 


/ 


which   shows    that    the    stream   lines   are   parabolic.      In  terms   of   the 
original   abscissa  x  =  ^a.  the   parabolas   are   given  by 


y^  -x.  ari  +  H^ 

U    TT 


TT     X        ,  TT        2         fr^„         ^\     X    '\ 

2  -?+2    T     -   (2    in   2)   - 

3> 


2         2-0- 
_   ,,   _  I  ,n22   +  - 


If  we  let  a  — ^  0  we  find  that  the  parabolas  reduce  to 


gx 

^1  =  -  7^ 

2u 


which  gives  the  shape  of  a  thin  heavy  sheet  of  fluid  projected 
horizontally  from  the  brink  with  velocity  u. 

The  first  order  approximation  to  the  equation  for  the  lower 
nappe  of  the  overfall  is 


(5.7) 


li^^Re 
a    TT 


£n 


i  1-  e~ 


ttA 


-  1 


Jl  -  e" 


dA  , 


0   <   i 


ttA 


+  1 


If  we  use  the  substitution  ^1  -  e  ^  =   o,    the  representation  (3.7) 
becomes 
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jt:p^ 


which  shows  that  the  equation  of  the  lower  nappe  can  be  expanded  in 
powers  of  yjl  -  e"""^^  if  0  ^  ^  <  oo  .   We  find  that  for  the  neighborhood 
of  the  brink  the  first  order  approximation  for  the  lower  nappe  is 


r^     ^  ^  fi  (2n+l)/2 

TT    n^  Vfer 


(3  8)       yi_   ^ev^Av^   1  \(l-e-^^) 
^^•^''       «.  -  -  2  2.^  2-^  2krT  2H+T 


From  this  approximation,  the  slope  of  the  lower  nappe  is  zero  at  the 
brink  but  its  curvature  is  infinite  there. 

The  first  order  approximation  for  the  brink  depth  ratio 

I  =   f(0,l;£) 
is 


(3.9) 


7^{0,1;b) 


This  is  the  same  as 
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yi(0,l;e) 


-i       V  Jl-e-^^  +  1 


dA 


1   +  -^  Re 

IT 


in 


7l-e-^^-l 


-oo 


Jl  -  e-^' 


dA 


+  1 


y 


1   +  ^ 


=     1    +    -^ 
TT 


0 

/ 

-oo 
0 


in 


Jl+  e-^^  -  1 


dA 


Jl+e-^Vl. 


ATrdA 


J^T^ 


-00     >iH-  e 
oo 


=  1  - 


~? 


in   ada 


IT        -,  a/l+a 


1 


37?  I 


CD     (-l)"r(n  +  ^)     r^ 


IT  '      n=0 


TTT 


1 


in  ada 
n+(3/2) 


1  - 


^ 


oo    (-l)'^r(n+-|) 


TT^/^  n=0     nl(n+-^) 


2 


(3.10) 


7^(0,  l;e) 


a 


=  1  -   e(.386) 


With  this,    the   first   order  equation  for   the   upper  nappe   is 
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y 


i+i 


—  =  1  +  ^  Re 
a      IT 


0 


.Jl-e-^^  +1, 


/"  i 


=  1+  -^  Re  / 

T  \ 


-ttA 


0     vJl-e-^'^+1. 


>/l-e-^^  -1 
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The  contraction  coefficient  b/a  for  the  case  e  =  1  is  important 
because  it  appears  in  various  formulas  for  flow  measurement  when  the 
overfall  is  used  as  a  metering  device  in  a  channel.   We  cannot  expect 
the  approximation  (3.I0)  for  the  brink  depth  ratio  to  be  good  if  e 
is  not  small.   However,  it  is  interesting  to  see  what  it  gives  if  the 
upstream  velocity  is  critical.   If  e  =  1,  (3«10)  gives 


yi(o,i;i) 


a 


=  .614  . 


This  is  to  be  compared  with  the  value 


(I)    =  .708 


which  is  the  average  of  contraction  coefficient  values 


19 


.715  Rouse  [2] 

.705  Southwell  and  Vaisey  [3] 

.720  Jaeger  [4] 

.720  Roy  [5] 

.710  Fraser  [6] 

.676  Hay  and  Marklemd  [7] 

deduced  from  various  experimental  and  mathematical  approximation  pro- 
cedures different  from  the  one  we  are  using  here.   Each  of  the  above 
papers  covers  the  case  of  critical  flow  only.   It  should  also  be 
remarked  that  some  of  the  values  in  the  above  table  are  not  approxima- 
tions to  the  actual  brink  depth  ratio  but  are  approximations  to  the 
stream  depth  ratio  corresponding  to  x  =  e,  where  e-u  is  small.   When 
this  is  the  case  the  corresponding  paper  does  not  contain  an  estimate 
of  just  how  small  e,  is. 

The  formula  (3. 10)  for  the  brink  depth  ratio  can  be  improved 
without  actually  finding  the  second  order  approximation  fp.   One  way 
to  accomplish  this  is  to  replace  f{0,T\;e.)   with 

(3.12)  f*(Ti;£)  =  c^+c^(r)-l)+-^  (ti-1)2 

where  the  coefficients  are  chosen  so  that 

f*(0;e)  =   f(0,0;e) 
f*(0;e)  =  f^(0,0;e) 

which  are  known  values;  and  so  that 
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(5.13) 


1+e'^f'^  (0,1) 
f*(l;e)  -  ' ^ 


1  -2e^f3_(0,l) 


The  right  member  of  (3. 13)  comes  from  the  boundary  condition  (2.13 ), 
namely. 


(3.14) 


2.  .    _         l+f^(0,l;e) 

\^^'^'^^   -   l+2e-  2ef(0,l;e) 


after  f(^,'rj;e)  in  the  right-hand  side  of  (3.14)  is  replaced  by 


Using 


i{i,r\;z)   =  T]  +efi(e,Ti) 


f,,(0,l)  =iRe  U^^^-^ 

./^T^^  +  1 


C=i 


and 


1  «   /^-  1 
—  Hn     i 


/2 +1 


f-^(0,l)  =  -.386 


a  computation  of  the  right  side  of  (3.13)  shows 


(l;,)  =  l^^e^('^l^)   . 

1+  £^(.772) 


The  value  f  (0, 0;e)  also  comes  from  (2.13)  and  it  is 


f  (0,0;£)  =  -^— 
"1         /I+2I 


The  coefficients  in  (3.12)  are  then  determined  by 
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f*(0;e)   =  c^-  c^+^  =   0  , 


f:;;(o;e) 


c-L-c^  = 


/l  +  2e 


f;(i;s) 


11+  e^(.3l4) 
l+e''^(.772) 


from  which   the    improved   approximation  to   the  brink  depth  ratio   is 


f*(i;8)  =  c^  =  c^  =^, 


'1    "    2 


/1+  2e, 


1 
2 


.vAT^ 


+  c. 


(3.15) 


f*(l;£) 


l  +  e 


.yi+2£       J  l+e 


'^(.314) 


''^(.772) 


This  result  can  be  expected  to  be  a  better  approximation,  for  a  wider 
ramge  of  values  of  s,  than  the  formula  (3 .10),  namely. 


yi(0,l;E) 


a 


=  1-  e(.386)  . 


For  e  =  .1  formula  (;5.10)  gives 


while  (3.15)  gives 


For  e  =  '5  we  have 


7^(0,  l;.l) 


=  .961 


f*(l;.l)  =  .955 
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YnfO,  1;.5) 

— =  .807 


f*(l;.5)  =  .829 


For  e  =  .9  which  means  that  the  flow  is  supercritical  to  the  extent 


that  u  =  '>i^—  /ga  =  1.05  yga  we  find 


y-l(0,l;.9) 
^^ =  -^53 

f*(l;.9)  =  .738  . 

For  e  =  1,  corresponding  to  critical  speed  we  have 

7.(0,1;!) 

^ -  .614 

a 

while 

f*(l;l)  =  .719  . 

This  last  result  compares  very  favorably  with  the  aver 

(|)  =    .708 

of  brink  depth  ratios  from  other  sources  cited  above. 

Better  approximations  can  presumably  be  found  by  passing  to  a 
determination  of  f^  but  for  the  time  being  the  results  presented 
above  appear  to  be  sufficiently  accurate. 
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